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In two space dimensions, the percolation point of the pure-site clusters of the Ising model coincides 
with the critical point Tc of the thermal transition and the percolation exponents belong to a special 
universality class. By introducing a bond probability pb < 1, the corresponding site-bond clusters 
keep on percolating at Tc and the exponents do not change, until = Pck = 1 — exp(— 2J/fcr): for 
this special expression of the bond weight the critical percolation exponents switch to the 2D Ising 
universality class. We show here that the result is valid for a wide class of bidimensional models 
with a continuous magnetization transition: there is a critical bond probability pc such that, for 
any pb >Pc, the onset of percolation of the site-bond clusters coincides with the critical point of the 
thermal transition. The percolation exponents are the same for Pc < Pa < 1 but, for pb = Pc, they 
suddenly change to the thermal exponents, so that the corresponding clusters are critical droplets 
of the phase transition. Our result is based on Monte Carlo simulations of various systems near 
criticality. 

PACS numbers: 64.60.Ak, 75.10.Hk 
Keywords: Phase transition, percolation 



I. INTRODUCTION 

The purely geometric percolation phenomenon I)] has 
astonishing analogies with ordinary second-order thermal 
phase transitions: the percolation variables show power- 
law behaviour near criticality, with relative exponents; 
the exponents are not all independent, but related to each 
other by scaling equalities which are formally identical 
to the corresponding thermal ones; the critical indices 
do not depend on the specific features of a system but 
they are grouped in universality classes, which in the 
percolation case are uniquely specified by the number of 
space dimensions of the system. 

These analogies suggested many years ago that a con- 
tinuous thermal transition might be nothing but a per- 
colation transition ^] , provided one defines suitable clus- 
ters, or droplets. The growth of the droplets describes 
naturally the propagation of correlations between differ- 
ent particles of the system and the formation of an in- 
finite spanning structure represents the long-range order 
of the system in the new phase. 

In this paper we focus on spin models with a continu- 
ous magnetization transition. In order to map the criti- 
cal behaviour of a system onto a percolation picture, one 
must set a correspondence between thermal and geomet- 
ric variables. The main percolation variables are: 

• the percolation strength P, i.e. the probability that 
a site chosen at random belongs to a percolating 
cluster; 



the average cluster size S, 



S 



(1) 



where is the number of clusters with s sites and 
the sums exclude eventual percolating clusters. 

The conditions for a successful mapping are listed below: 



• the percolation point must coincide with the ther- 
mal critical point; 

• the connectedness length (average cluster radius) 
diverges as the thermal correlation length (same 
exponent); 

• the percolation strength P near the threshold varies 
like the order parameter m of the model (same ex- 
ponent); 

• the average cluster size S diverges as the physical 
susceptibility x (same exponent). 

It is not yet known whether, given a system, there are 
at all clusters satisfying such conditions, which are quite 
strict. The first studies concentrated on the simplest the- 
ory, the Ising model without field. A real breakthrough 
in these investigations was the discovery that the Ising 
model can be rewritten as a geometrical model 0] , where 
the fundamental objects are special site-bond clusters, 
i.e. clusters built by joining nearest-neighbouring spins of 
the same sign with a temperature-dependent bond prob- 
ability Pck = 1 — exp(— 2J/fcr) (J is the Ising spin-spin 
coupling). It is possible to prove that these site-bond 
clusters are just the critical droplets of the system. Start- 
ing from this result, which is valid more in general for the 
q-state Potts model, it is possible to define the droplets 
for a wide variety of systems, like models with several 
spin-spin couplings, as long as they are all ferromagnetic 
Q, and 0{n) spin models A common feature of all 
these results is the fact that a geometrical bond (with a 
bond probability) is associated with each spin-spin cou- 
pling in the percolation picture. In this way, for mod- 
els with several (ferromagnetic) interactions, the clusters 
are quite weird objects, consisting of geometrically dis- 
connected parts which are joined to each other by virtue 
of invisible bridges due to long-ranged interactions. Be- 
sides, if not all couplings are ferromagnetic, it is not yet 
clear if and how one can define the droplets. 
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Very recent results show that the simple connections 
between nearest-neighbour spins have a close relation- 
ship with the phase transition of several bidimensional 
models Q: the percolation point of the pure-site clus- 
ters coincides with the critical point and the perco- 
lation exponents, which differ from the thermal expo- 
nents, are the same for models in the same universal- 
ity class. In our opinion this property indicates that, in 
many cases, the simple geometrical connectivity between 
nearest-neighbours plays a crucial role in the mechanism 
of the phase transition. 

In the 2D Ising model, the fact that both the pure-site 
and the site-bond clusters with bond probability pcK 
start to percolate at the critical temperature Tc neces- 
sarily implies that, by taking a bond probability ps such 
that PCK < Pb < 1, the percolation temperature of the 
site-bond clusters with bond weight is always Tc. A 
renormalization group analysis led to the conclusion 
that the percolation exponents are the same as for the 
pure-site clusters for any ps > PcK , and that they sud- 
denly switch to the thermal exponents for pb = Pck- 
For Pb < Pck, the site-bond clusters percolate at some 
temperature Tp < Tc and the exponents switch to the 
2D random percolation universality class. This analysis, 
though reliable, is however not rigorous. In this paper we 
will show that the result is true and valid for many bidi- 
mensional models. We will see that, given a model, there 
is a critical probability pc such that, by taking a bond 
weight Pb > Pc (including the pure-site case pb — ^) the 
site-bond clusters keep percolating at Tc and the critical 
percolation exponents are the same; for = pc the ex- 
ponents change abruptly to the thermal exponents. The 
probability Pc, which is model-dependent, is the minimal 
probability which still allows the formation of a percolat- 
ing cluster at Tc, and this gives us a general prescription 
for the critical droplets. 



II. RESULTS OF THE SIMULATIONS 

Our investigations consisted of Monte Carlo simula- 
tions on square lattices of several models around the crit- 
ical temperature. We took the same models that were 
studied in i.e. 

1. the Ising model, 7i = —JJ^ij^i^j > 0, Si — 

±1); 

2. a model with nearest-neighbour (NN) ferromag- 
netic coupling and a weaker next-to-nearest (NTN) 
antiferromagnetic coupling: Ti. = — Ji J2nn ^i^j ~ 
J-iY^NTN^^Sj (-^1 > 0, J2 < 0, IJ2/J1I = 1/10, 
Si = ±l); 

3. the continuous Ising model, Ti = -~JJ2ij^i^j 
(J>0, -l<5i< + l); 

4. SU(2) pure gauge theory in 2+1 dimensions; 



5. the 3-state Potts model, H = —JJ2ijH^i,Sj) 
(J>0, = 1,2,3); 

6. a model obtained by adding to 5) a weaker next- 
to-nearest (NTN) antiferromagnetic coupling: H = 

J2 < 0, IJ2/J1I = 1/10, s, = 1,2,3); 

The critical temperatures of all these systems are known 
with good precision (see |^), for Models 1 and 5 the val- 
ues are analytically known. We stress that the models 
belong to the 2D Z{2) (1-4) and Z(3) (5,6) universality 
classes. The estimates of Q for the main critical indices 
of the percolation transition of the pure-site clusters for 
the two classes are in very good agreement with theoret- 
ical predictions ^st these analytical values in 
Table n 
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153/80 


0.932(2) 



TABLE I: Critical percolation indices for pure-site clusters 
for the two universality classes we considered. 



For the simulations we made use of standard algo- 
rithms like Metropolis or heat bath; in some cases we 
could apply the Wolff cluster algorithm, which allowed 
us to reduce sensibly the correlation of the data. For 
each system we took four to six different lattice sizes. 
The clusters were identified by means of the algorithm 
devised by Hoshen and Kopelman We used every- 
where free boundary conditions for the cluster labeling 
and say that a cluster percolates if it connects the top 
with the bottom side of the lattice. At each iteration we 
calculated the variables P and S and the size Sm of the 
largest cluster, from which one can determine the fractal 
dimension D of the percolating cluster at Tc- Moreover, 
from the data sample of the percolation strength P one 
can extract a variable which turns out to be a formidable 
tool to investigate numerically the percolation transition. 
In fact, the number of configurations with (at least) a 
percolating cluster (i.e. for which Pj^Q) for a given tem- 
perature and lattice size, divided by the total number of 
configurations, returns a scaling variable 11 called perco- 
lation cumulant, whose properties are identical to those 
of the Binder cumulant in standard thermal transitions 
[T^ . In particular, the value of 11 at the critical threshold 
is a universal quantity, i.e. it is the same for models in the 
same universality class. We determined the percolation 
exponents through standard finite size scaling tecniques 
at the critical temperature Tp, considering simply the 
leading behaviour 
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P{Tp) oc L-^"/"" (2) 
S{Tp) oc L'^"/'"'' (3) 
Sm{Tp) (X i^. (4) 

We began our studies starting from the Ising model. 
We tested three different expressions for the bond prob- 
abihty [3 Pb- CKi = 1 - cxp(-3 J/ZcT), CK2 = 
1 - exp(-2. 5 J/fcT), Ci^a = 1 - exp(-2.2 J/fcT). In all 
cases we found that, even for the smallest lattice size we 
examined (100^), the threshold value of the percolation 
cumulant 11 did not differ from the value relative to the 
pure-site clusters (0.9832(4)), which is a clear indication 
that the critical indices of the three percolation tran- 
sitions remain in the universality class of the 2D Z{2) 
pure-site clusters. Figs. ^ and |21 show the situation for 
the probabilities CK2 and CK^, respectively. 
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FIG. 1: 2D Ising model: value of the percolation cumulant 11 
at the critical temperature Tc for various lattices. The bond 
weight is CK2 = 1 — exp(— 2.5J/fcT). All values agree with 
each other (so Tp — Tc) and with the (universal) pure-site 
value, represented with its error by the dashed lines in the 
plot. 

Indeed, the values of the percolation exponents we ex- 
tracted from finite size scaling fits of P, S and Sm at the 
critical point are in good agreement with the Z{2) val- 
ues reported in Table HI Next, we approached the critical 
value of pcK at Tc {pck{Tc) — 0.58578...) by successive 
numerical trials. We found that, up to a value greater 
than pck{Tc) by 4 — 5%, the threshold value of the cumu- 
lant did not change for any of the lattices we have taken, 
and the situation is the same as in Figs. ^ and |21 For 
still lower probabilities the cumulant decreases apprecia- 
bly for the smaller lattices and the less the larger the size. 
That is the " lattice reaction" to the fact that we are near 
an abrupt discontinuity. We know that the finite size of 
the lattice reduces divergences to finite peaks and discon- 
tinuities to smooth variations. Close to a sharp disconti- 
nuity, the real behaviour of the system is reproduced only 
on very large lattices; for small sizes, the system feels the 
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FIG. 2; 2D Ising model: value of the percolation cumulant 11 
at the critical temperature Tc for various lattices. The bond 
weight is CK3 = 1 — exp(— 2.2 J/fcT). All values agree with 
each other (so Tp — Tc) and with the (universal) pure-site 
value, represented with its error by the dashed lines in the 
plot. 



proximity of the new state and that gives rise to interme- 
diate configurations, which are unphysical. On the other 
hand, the threshold value of the percolation cumulant on 
our largest lattice (1000^) remained fixed at the pure-site 
value until pb/pck{Tc) « 1.02, and that clearly suggests 
that the transition from the pure-site behaviour to the 
droplet behaviour takes place just when pB = Pck{Tc)- 

By taking exactly ps — Pck{Tc), the threshold value 
of the percolation cumulant is again the same for all lat- 
tices, and equal to the value 0.585(1) that labels the 
universality class of the 2D Z{2) droplets. Next, we 
took some values of pb < Pck{Tc)- Just below pck{Tc) 
we expected to see the same mixed behaviour observed 
when one approaches the critical probability from above, 
and that is indeed the case. However, already for 
Pb/pck{Tc) < 0.96 we could see again a "clean" be- 
haviour of the percolation variables on all lattices: the 
percolation cumulants no longer cross at Tc but at a 
somewhat lower temperature Tp, which decreases by de- 
creasing Pb ■ The value of the cumulant at the crossing 
point is the same for any pB and agrees with the value 
corresponding to the 2D random percolation universality 
class (0.450(2)), as shown in Fig. 13 Our findings confirm 
the scenario predicted in 7]. 

As far as the other models with an Ising-like transition 
are concerned we remark that, except for the continuous 
Ising model jM] , the problem of the existence of eventual 
critical droplets is open. Since in each case the pure-site 
clusters begin to percolate at Tc, it was natural for us to 
check whether also here one has the same situation like 
in the Ising model. We then sat at Tc in every case, and 
introduced progressively decreasing bond weights. We 
remark that for the continuous spin models 3 and 4 we 
bound nearest-neighbour spins of the same sign, indepen- 
dently of their absolute value. 
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FIG. 3: 2D Ising model: value of 11 at the percolation 
temperature Tp for various lattices. The bond weight is 
1 — exp{—J/kT) < pc = 1 — exp(— 2J/fcT). In this case is 
Tp < Tc- The data points agree with the (universal) value 
for 2D random percolation, represented with its error by the 
dashed lines in the plot. 
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2D Ising 
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0.585(1) 


Model 2 


0.131(10) 


1.742(12) 


1.862(20) 


0.583(4) 


Model 3 


0.121(9) 


1.764(14) 


1.870(11) 


0.587(3) 


Model 4 


0.140(19) 


1.761(17) 


1.882(18) 


0.586(5) 



TABLE II: Critical percolation indices for the site-bond clus- 
ters of the models with an Ising-like transition when pb = Pc, 
compared with the values of the 2D Ising droplets. 



We found that the scenario is indeed the same. In ev- 
ery case, we found a critical probability pc, which equals 
0.583(1), 0.6115(9) and 0.6275(7) for Model 2,3 and 4, re- 
spectively. We adopted several values for ps, both above 
and below pc- For > pc we have always found that the 
percolation point coincides with the critical temperature 
Tc of the model, and that the critical percolation indices 
coincide (within errors) with the Z(2) percolation indices 
of Table m In all cases we could show analogous pictures 
as Fig. n or 121 If < Pc, the percolation temperature 
Tp < Tc and the exponents belong to the 2D random per- 
colation universality class (as in Fig. O. Like in the Ising 
model Pc is then the minimal probability for which the 
site-bond clusters still percolate at Tc. Fig. ^ shows the 
threshold value of the percolation cumulant for all Z{2) 
models when ps — Pc] all values agree within errors. In 
this case, the critical indices are in accord with the ones 
of the 2D Ising droplets, as shown in Table UTI 

Finally we investigated Models 5 and 6. We know that 
the g-state Potts model can be transformed into a cluster 
model by means of the Fortuin-Kasteleyn transformation 
0. The clusters are defined in the same way as for Ising, 



FIG. 4: Threshold value of the percolation cumulant 11 at the 
critical temperature Tc by using the "minimal" probability pc 
for all Z{2) models. The X-coordinates 1, 2, 3 and 4 label the 
four models we investigated. There is a clear accord of the 
four data points with each other. The dashed lines indicate 
the bounds of the more precisely measured Ising value (first 
point to the left), which labels the universality class of the 2D 
Ising droplets. 



except that the bond weight pfk = 1 — exp(— J/fcT), but 
the situation is identical. It is then not surprising that for 
the 2D (/-state Potts model the same scenario as for 
Ising was predicted |lO|. The results of our simulations 
confirm such prediction. The situation does not change 
for Model 6: we found a critical bond probability pc — 
0.61(1), at which the site-bond clusters become critical 
droplets for the model, as one can see in Table IIIII In 
Fig. [SI we compare the threshold values of 11 for the two 
Z(3) models when pB = Pe- 
rn 0.68 I ^ ^ ^ . . 1 
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FIG. 5: Threshold value of the percolation cumulant 11 at the 
critical temperature Tc by using the "minimal" probability pc 
for the Z{3) models. The X-coordinates 5 and 6 label the two 
models we investigated. The values agree well within errors. 
The dashed lines indicate the bounds of the more precisely 
measured 3-state Potts value (point to the left), which labels 
the universality class of the 2D Z{3) droplets. 
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26/15 


28/15 


0.649(9) 
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0.143(17) 


1.725(21) 


1.858(18) 
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TABLE III: Critical percolation indices for the site-bond clus- 
ters of Model 6 when p_B = Pc, compared with the correspond- 
ing values of the 2D 3-state Potts droplets. 



III. CONCLUSIONS 

We found a general criterion to identify the critical 
droplets of various bidimensional systems. The droplets 
are the clusters obtained by joining nearest-neighbour 
spins of the same sign (-state for Potts-like systems) 
with the minimal bond probability Pc that still allows 
the clusters to percolate at the critical point of the ther- 
mal transition. For any ps > Pc the onset of percola- 
tion for the relative site-bond clusters keeps coinciding 
with the thermal threshold and the critical indices are 
the same up to the pure-site case pB = I- We stress 
that this prescription seems to work also for the cases 
in which one can rigorously define the droplets. For the 
continuous Ising model, for instance, the droplets are de- 
fined by introducing local bond weights which also de- 
pend on the length of the spins and not only on their 
sign These clusters certainly differ in detail from 

the droplets defined in this paper; nevertheless their be- 
haviour at criticality is identical. The matter would be 



even more involved for models with several ferromagnetic 
interactions, for which a rigorous definition of droplets re- 
quires the presence of longer-range connections than just 
between nearest-neighbour spins The fact that the 
critical behaviour of the system is reproduced by simple 
site-bond clusters, no matter how complicated the model 
is, suggests that the long-range fluctuations of the sys- 
tem, which are responsible for the phase transition, are 
embodied in the simple magnetic domains of the system; 
the bond probability Pc is necessary in order to destroy 
additional spin correlations due to purely geometrical ef- 
fects, as already remarked in Besides, we have seen 
that the result is also valid in cases where a rigorous def- 
inition of droplets is, at present, missing (Models 2, 4 
and 6). In particular, it remains true for SU{2) pure 
gauge theory, which is a very complicated model involv- 
ing many different interactions, like multispin couplings, 
long-ranged couplings (ferromagnetic and antiferromag- 
netic) and self-interactions. This shows that the result 
has some generality, and it would be interesting to check 
to which extent it is true for bidimensional models with a 
continuous magnetization transition. It would be also in- 
teresting to check whether simple site-bond clusters play 
a role in the description of critical behaviour at higher 
dimensions as well. 
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